Thermodynamics and phase transitions for the Heisenberg model on the pinwheel 

distorted kagome lattice 
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We study the Heisenberg model on the pinwheel distorted kagome lattice as observed in the 
material Rb2Cu3SnFi2. Experimentally relevant thermodynamic properties at finite temperatures 
are computed utilizing numerical linked-cluster expansions. We also develop a Lanczos-based zero- 
temperature numerical linked-cluster expansion to study the approach of the pinwheel distorted lat- 
tice to the uniform kagome-lattice Heisenberg model. We find strong evidence for a phase transition 
before the uniform limit is reached, implying that the ground state of the kagome-lattice Heisenberg 
model is likely not pinwheel dimerized and is stable to finite pinwheel-dimerizing perturbations. 
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I. INTRODUCTION 

Among all of the realistic quantum spin models studied 
to date, the spin-half kagome-lattice Heisenberg model 
(KLHM) stands out as one with perhaps the most es- 
oteric low-energy properties^ The classical Heisenberg 
model on this lattice of corner-sharing triangles has an 
exponentially large number of ground states. Fluctua- 
tions lead to a selection of \/?> x order.— However, 
such an order is clearly absent in the spin-half models A 
comparison between exact diagonalization (ED) results 
in different lattice geometries highlights the unusual spin 
physics of this modeli Square and triangular lattice mod- 
els show a tower of rotor states separated from the rest 
of the spectrum, as expected for a system with sponta- 
neously broken spin-rotational symmetry. In contrast, 
the KLHM shows that (i) the energy gaps are an or- 
der of magnitude smaller than on the other lattices, and 
(ii) there are a very large number of low-lying singlet 
states. These data have been interpreted to support the 
existence of (i) a very small or even zero triplet gapj^ 
(ii) gapless singlet excitations^ and (iii) an exponentially 
large number of singlets below the lowest triplet^ 

Numerical and analytical studies have suggested al- 
ternative ground states, which include several different 
valence-bond crystals (VBCs), as well as gapped and 
gapless resonating valence-bond (RVB) spin liquids**"— 
Most numerical studies, including the recent density ma- 
trix renormalization group (DMRG) study of Yan et 
al-r^. suggest a singlet ground state with a gap to spin 
excitations in the thermodynamic limit. While earlier 
numerical studies^ provided evidence for a spin gap of 
the order of J/20, where J is the nearest-neighbor ex- 
change constant, and a much smaller gap to singlet exci- 
tations (possibly even gapless) , recent DMRG work puts 
the singlet gap at 0.05 J, with a larger gap for the triplets. 
These studies should be contrasted with the experimental 
results for the KLHMJ^r— Among them, the herbert- 
smithites with structurally perfect kagome planes have 
attracted the most interest. The experimental data for 
these systems show no evidence of a spin gap. How- 



ever, substitutional impurities and Dzyaloshinskii-Moria 
anisotropy are present in these systems and may be 
playing a role in producing the eventual thermodynamic 
phasej 14 ' 15 

One kagome material that does show gapped spin exci- 
tations is Rb2Cu3SnFi2i 16 ' 17 This material also has sub- 
stantial pinwheel distortions, which results in a 12-sitc 
unit cell^ and four different nearest-neighbor exchange 
constants, with the ratio of the largest to smallest being 
about two (see Fig. [IJ. Given the many theoretical pro- 
posals for spin-gapped and VBC phases for the KLHM 
and this experimental finding, it is natural to ask if the 
uniform KLHM has a tendency to spontaneously dimer- 
ize in this pattern. 

Here, we study this possibility and the role such a pat- 
tern of distortion might play in the ground state of the 
uniform KLHM. In earlier series-expansion studies, dif- 
ferent dimer configurations were studied order by order in 
powers of a parameter A, which takes one from a dimer- 




FIG. 1. (Color online) Kagome lattice with bond anisotropy 
as relevant to the Rb2Cu3SnFi2 material. The strength of 
the bonds (J a ) can be inferred by their thickness. Red (thick) 
bonds show dimers of the pinwheel VBC order with a 12-site 
unit cell, which is the ground state when J\ is sufficiently 
larger than the other exchange interactions. The dotted- 
dashed lines show the 12-site and 24-site periodic clusters used 
in ED. 
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ized to a uniform Hamiltonian4 In each successive or- 
der, the set of dimer states was reduced by keeping only 
those that were lowest in energy. In low powers of A, the 
pinwheel-dimer state has one of the highest energies, as 
it has no resonating hexagons, and was not considered 
further.— However, all dimer states were found to be de- 
generate to the second order in perturbation theory and 
stayed within a few percent of each other for A = 1 to 
the highest order studied. Thus, it is possible that higher 
powers of A can stabilize a different state than the ones 
picked out perturbatively. Standard series expansions in 
A are clearly not suitable to study such a selection. 

We develop here a numerical linkcd-clustcr expansion 
(NLCE}i£ to study this problem. This is a graphical 
method based on the principle of inclusion and exclusion 
for calculating the properties of an infinite system by 
summing over contributions from all smaller clusters. Its 
convergence is not based on the existence of a small pa- 
rameter, but rather on having a short correlation length, 
so that contributions from larger clusters become small. 
Such an expansion always converges at high tempera- 
tures, but can also converge at low temperatures for mod- 
els where the correlation length remains shorts 

A suitably chosen NLCE scheme can also converge at 
T = if the system has long-range VBC order, but the 
fluctuating correlation length on top of the VBC order 
remains small. By building the clusters in terms of the 
dimers of the pinwheel VBC order (see Appendix for de- 
tails of the expansion), the properties of such a phase can 
be calculated. There- is no small parameter in the expan- 
sion, as the properties of the finite clusters are calculated 
exactly. For concreteness, we introduce a family of mod- 
els that depend on a parameter A, defined in the same 
way as for the series expansion. As long as our NLCE 
shows internal convergence, one can study any ground- 
state property. The failure to converge signals a phase 
transition away from the VBC phase. 

We first use the NLCE to study the specific heat (C v ), 
entropy (S), and spin susceptibility (x) of the pinwheel 
distorted kagome lattice with exchange constants that de- 
scribe Rb2Cu3SnFi2 - 16,17 Our results are accurate in the 
thermodynamic limit and converge at almost all temper- 
atures (including T = 0). They also compare very well 
with the experimental data for the uniform susceptibil- 
ity, confirming the pinwheel valence-bond nature of the 
ground state of this material. We then consider a model 
in which the exchange constants are unity for the strong 
bonds in the pinwheel order and A for all others. We 
study the thermodynamics of this model for A from the 
dimerized (A <C 1) to the uniform KLHM (A = 1). Since, 
sufficiently away from A = 1, the NLCE converges at 
T = 0, we use a Lanczos-based NLCE and include larger 
cluster sizes to study ground-state properties (see Ap- 
pendix for details) . We find that a VBC order parameter 
decreases as A increases and, according to a power-law 
fit, vanishes at A c < 1. The results for finite-size clusters 
with periodic boundary conditions point to a first-order 
phase transition from the VBC phase to a spin-liquid-likc 




FIG. 2. (Color online) Specific heat, entropy, and uniform 
spin susceptibility per site for the pinwheel distorted kagome 
lattice of Fig. [TJ with (a)-(c) J 2 = 0.95 Ji, J 3 = 0.85Ji, and 
J 4 = 0.55Ji (the model for Rb 2 Cu 3 SnFi 2 ), and (d)-(l) J 2 = 
J[i = Ji = AJi. The last two orders of the bare sums in the 
NLCE are shown. They are 5 th and 6 th orders for A = 0.50, 
and 8 th and 9 th orders for A = 0.75 and A = 1.00, and in 
(a)-(c). Dashed lines are results from ED for the periodic 12- 
site cluster. Blue circles in (c) are experimental results from 
Ref.QJ (see text). 

phase at A c > 0.9. 

II. THE MODEL 

The KLHM Hamiltonian is written as: 

H = ]T J a s, • S, (1) 

where (..) denotes nearest neighbors, Si is the spin-1/2 
vector on site i, and a is the bond type. We consider the 
bond anisotropy seen in Fig. Q] with four distinct cv's. 

III. RESULTS 

A. Thermodynamics of Rb 2 CuaSnFi 2 

We study the thermodynamics of the anisotropic 
KLHM shown in Fig. [TJ with the exchange constants sug- 
gested for Rb2Cu3SnFi2 in a recent experiment*^ We 
take J2 = 0.95Ji, J3 = 0.85Ji, and J4 = 0.55Ji, and 
set the unit of energy to J\. In Figs. [Ha)-[2jc), we show 
results for the 8 th and 9 th orders of bare NLCE sums^ 
along with results from ED on a 12-site periodic clus- 
ter, for C v , S, and \ per site. In addition to the high- 
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T region, the thermodynamic quantities also fully con- 
verge in a low-T window between and 0.1. This shows 
that clusters in our dimcr expansion around the pinwheel 
VBC order are able to capture the low-T properties very 
well, which in turn, supports the finding that this ordered 
phase describes the ground state of Rb2Cu3SnFi2.— ^ 

In Fig.[2ja), note the relatively flat C v of this model for 
T between 0.1 and 1.0 [Fig. dta)], reminiscent of a classi- 
cal gas. The entropy in that region varies almost linearly 
with logT [Fig. EJb)] . The ED results for C v exhibit 
a double-peak structure but clearly underestimate the 
value of the high-T peak as seen in Fig.[^a). Large finite- 
size effects in ED can also be seen at low T in Fig. EJc). 

Figure [2jc) shows a comparison between the exper- 
imentally measured^ x an d the NLCE results using 
Ji = 234 K and the g factor, g = 2.46 (needed to 
match the two at high T). We find a remarkable agree- 
ment between the NLCE results within its convergence 
region (T > 0.4 and < T < 0.15) and the experi- 
ments [Fig. EtcYL The fact that the zcro-T fit of the 
spectra in Ref. [l7| and the finite- T susceptibility corre- 
spond to roughly the same parameters also shows that 
the exchange constants are nearly temperature indepen- 
dent. Small deviations at T < 0.07 can be attributed to 
the uncertainty in the impurity subtraction^ 



B. Thermodynamics of the model with 

J% = Jz = Ja = AJi 

We are interested in the behavior of thermodynamic 
quantities and the NLCE convergence at low T for dif- 
ferent A valuesJ^ In Figs. HKd)-|Hl), we show C v , S, and 
X for A = 0.50, 0.75, and 1.00. Figures |3(d)-^f) show 
excellent convergence of the NLCE results at all T in the 
5 th and 6 th orders for A = 0.5. Note that the specific 
heat in Fig. HJd) has a single peak at T ~ 0.27 with no 
other features, which is characteristic of a single dimcr. 

Upon increasing A to 0.75, the NLCE results still con- 
verge at nearly all temperatures by increasing the order 
to 9. Despite this convergence, which suggests that the 
physics of the VBC is still dominant, the temperature 
dependence of the specific heat, entropy, and suscepti- 
bility changes from those of A = 0.5. The peak in the 
specific heat splits into two peaks [Fig. [Sfg)] an d the en- 
tropy develops a hump in that region [Fig. HJh)]. Also, 
the peak in susceptibility broadens and moves to lower 
T [Fig. [210] ■ Unlike the A = 0.5 case, due to finite-size 
effects, the 12-site ED results do not completely agree 
with the NLCE results when A = 0.75, and again, under- 
estimate the specific heat for 0.2 < T < 0.5. 

The convergence of the NLCE results changes com- 
pletely in the uniform limit (A = 1), as depicted in 
Figs. [2Jj)-[2Jl). Similar to previous NLCE calculations 
that use the site or triangular expansions on the kagome 
lattice^ the series do not converge below T ~ 0.3. 
Hence, the low-temperature physics can no longer be de- 
scribed in the dimcrized picture. 
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FIG. 3. (Color online) VBC order parameter at T — vs the 
anisotropy parameter A for J2 = J3 = J4 = AJi. The re- 
sults for finite-size clusters show a first-order transition. Bare 
NLCE sums do not converge beyond A ~ 0.85. We use the 
Wynn sum with four cycles of improvement to extend the 
convergence to A ~ 0.9. Results for the last two orders of the 
former are fit to ^4(A C — A) s (thin dotted and dashed lines). 
Vertical solid lines show the region used for the fit. The inset 
shows the four bond energies for the 24-site ED vs A. 

C. The ground state 

To study ground-state properties directly, we use the 
Lanczos algorithm and carry out the NLCEs to the 13th 
order. We calculate the bond energies B a per bond of 
each type in Fig. [TJ The difference between B\ and any 
other bond energy, e.g., B2 — B±, can be used as the 
order parameter. For small A, when the system is al- 
most completely dimerized in the pinwheel pattern, B\ 
is dominant and the order parameter is maximum. As A 
increases, all B a should approach the same value, which 
at A = 1 is half the ground-state energy per site. Al- 
though this order parameter can never be exactly zero 
for any A ^ 1, there is the possibility of a phase transi- 
tion from the bond-solid order to some other phase for 
A c < 1, beyond which the difference in bond energies 
arises from the difference in the exchange constants, not 
dimerization. Naively, the order parameter would vanish 
linearly with A in this region. 

In Fig. [3J we show the dependence of B 2 — B\ on A. 
We find very good convergence for small values of A. 
However, for A > 0.85, the bare sums in the NLCE do 
not converge. Re-summing the series using the Wynn 
algorithm^ extends the region of convergence to A ~ 0.91 
(last two orders are shown in Fig. [3]). This is, in fact, an 
indication that a small pinwheel distortion of the kagome 
lattice fails to stabilize the VBC order and that there is a 
transition at 0.9 < A c < 1 to another phase [presumably 
a spin-liquid (SL) phase] in the thermodynamic limit. In 
Fig. [3l we also show the ED results from the 12- and 
24-site clusters for the entire region of A. Both clusters 
exhibit a clear first-order transition with A c < 1. In the 
inset of Fig. [31 we show the behavior of individual bond 
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FIG. 4. (Color online) Upper panel: Plot of the VBC order 
parameter B2 — £>i in the T-X plane for the 12-site ED cal- 
culations. The decrease in B2 — B\ by increasing A from the 
left VBC region to the right region is less abrupt at higher 
temperatures. Lower panel: The same quantity vs T for five 
different values of A below and above the transition point. 
Thin lines are the 7 th and 8 th order NLCE, and thick lines 
are the results for the 12-site cluster. 



energies for the 24-site cluster. It is interesting to see that 
the bond energy B3 of the so-called "empty triangles' — 
is the least affected by the transition. We should stress 
that B2 — B\ in the ED results exhibit large finite-size 
effects in the region where the NLCE is converged, and 
increasing the cluster size from 12 to 24 does not improve 
the value of B2 — B\ for A < A c . As one can see in Fig. [31 
for A < A c , the results from the 24-site cluster are fur- 
ther away from the (thermodynamic limit) NLCE results 
than those from the 12-site cluster. 

A power-law fit of the high-A NLCE results for the 
order parameter allows us to estimate A c in the event 
of a second-order phase transition. We find A c ~ 0.98 
when considering the last two orders of the NLCE after 
rcsummation. Results for the fit are shown as thin dotted 
and dashed lines in Fig. [3] If the transition is actually 
first order, A c will be smaller. 



D. Order parameter at finite temperatures 

Finally, in connection with the experiments, and to ex- 
plore the effects of a possible quantum critical region, we 
use NLCEs and ED on the 12-site cluster to examine the 
behavior of B2 — B\ at finite T. The contour plot in the 
upper panel of Fig. U shows the order parameter in the 
T — A plane obtained from the ED of the 12-site cluster. 
The blue area (dark area in the left) for A < 0.9 at low 
temperatures indicates the VBC phase, and the red area 
(dark area in the right) for A > 0.9 belongs to the other 
phase. Separating these two phases is a crossover re- 
gion, shown by the light color in that figure, which starts 
at the finite-size transition point of the 12-site cluster 
(A c ~ 0.91) and is broadened as it moves to smaller A 
at higher T. We expect the same qualitative behavior in 
the thermodynamic limit, with a transition point that is 
closer to the A = 1 axis. 

To compare with the NLCE results, in the lower panel 
of Fig. |H we plot the order parameter for different values 
of A below and above the transition point. An interest- 
ing feature there is the initial low-temperature rise of the 
order parameter above but very close to A c as T is in- 
creased from zero. In this region, thermal fluctuations 
initially increase the value of the order parameter before 
suppressing it at higher T. If the transition in the ther- 
modynamic limit is second order or weakly first order 
and closer to A = 1, there could be a quantum critical 
region that controls the temperature dependence of the 
KLHM over some temperature ranged 



IV. SUMMARY 

In summary, we have studied the KLHM on a pin- 
wheel distorted kagomc lattice using NLCEs and exact 
diagonalization of finite clusters. Wc have examined the 
thermodynamic properties of the model with exchange 
parameters appropriate for Rb2Cu3SnFi2, and obtained 
results that are consistent with the experimental obser- 
vations. We have also developed a Lanczos-based NLCE 
and studied the approach to the uniform KLHM at T = 
in such a pinwheel-dimerized system. We find strong ev- 
idence for a phase transition before the uniform limit 
is reached, implying that the KLHM may not sponta- 
neously dimerize in this pattern and its ground-state 
phase is stable to small pinwheel distortions. 
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APPENDIX: NUMERICAL LINKED-CLUSTER 
EXPANSION FOR THE DISTORTED KAGOME 
LATTICE 

In standard series expansions for the distorted kagomc- 
lattice Heisenberg models^ the ground-state properties 
of each cluster in the series are calculated using perturba- 
tion theory in the small parameter, A. The contributions 
of all of the clusters, up to a certain size, to each power 
of A are then added properly into the series. In contrast, 
in the NLCEs, the property of each cluster is computed 
exactly, i.e., to all orders of A, using exact diagonalization 
techniques J£ Because of this, the NLCE does not require 
A to be small, and one can achieve convergence in a region 
that extends beyond that of the usual series expansions. 
The convergence in NLCEs is controlled by the relevant 
correlation length in the system. Hence, unlike for stan- 
dard series expansions, by adding larger cluster sizes, we 
can always improve the convergence. 

In our expansion, the building blocks for generating 
the linked clusters are dimers of the pinwheel VBC. This 
means that in the nth order, we generate all the clusters 
with n dimers that can be embedded in the anisotropic 
kagome lattice (Fig. [TJ. Starting from a single dirtier 
in the first order, there are many possibilities for the 
topology of the clusters in the second order, generated 
by adding another neighboring dimer. However, among 
those, only the four clusters shown in Fig. [5] are topolog- 
ically distinct^. Next, we require that dimers be con- 
nected only through complete triangles. For example, in 
the second order, we consider only cluster D among those 
in Fig. [5j The above condition significantly reduces the 
number of clusters we have to diagonalize in higher orders 
while still capturing the fundamental physics at low ener- 
gies in relatively short expansions. This selection of the 
clusters offers a consistent scheme for the linked-clustcr 
expansion of the infinite lattice. Following the same pro- 
cess, we end up with two topologically distinct clusters 
in the third order by adding the next dimer, and so on. 
The clusters in the series, up to the fourth order, are de- 
picted in Fig. |H1 In Table. HI we have listed the number 
of topological clusters in each order up to the 13 th order. 

The expansion described above offers a unique ap- 
proach in which the series capture not only the high- 



A) B) C) D) 

_ / _\ _\ _\ 

J 1 J 3 J 2 

FIG. 5. (Color online) All possible configurations for two 
neighboring dimers (red thick bonds) on a kagome lattice with 
pinwheel anisotropy. In general, four different types of bonds, 
Ji, . . . , Ji, can be distinguished in this lattice. Cluster D, in 
which the two dimers are connected via a complete triangle, 
is the only cluster considered in the second order of our ex- 
pansion. 
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FIG. 6. (Color online) All the topologically distinct clusters 
up to the fourth order in our expansion. Bond types are the 
same as in Fig. [5] 



temperature physics, like the standard NLCEs used pre- 
viously , 18 ' 21 but also capture the low-temperature physics 
very well, as long as the pinwheel VBC is the ground 
state. This makes the NLCE used in this work fundamen- 
tally different from the ones implemented before. It uses 
a basis tailored to describe the ground state of this partic- 
ular model. Hence, we use it to calculate properties both 
at finite temperatures and, specially, at zero temperature. 
For finite-temperature calculations, we compute the full 
spectrum of the Hamiltonian for every cluster in the se- 
ries using standard full diagonalization techniques, and 
for zero-temperature calculations, we employ the Lanczos 
algorithm to access its ground state only. The advantage 
of the latter is that it enables us to deal with systems 
with much larger Hilbert spaces, and so we are able to 
carry out the expansion to significantly larger orders in 
comparison to the finite-temperature case. 

Computationally, the dominant part of the calculations 
is the diagonalization step in the last order. Ultimately, 
wc arc limited cither by the large number of clusters that 
exist in the largest order we consider and the compu- 
tational time it takes to solve for each of them, or by 
memory requirements for diagonalizing the Hamiltonian 
of each of those clusters. For the full diagonalization 
finite-temperature calculations, we consider clusters up 
to the 9th order, and for zero-temperature Lanczos cal- 
culations, up to the 13th order. Hence, only in the last 
order of the former case, 102 clusters need to be solved. 
Each of them has 18 sites and, therefore, the size of the 
Hilbert space for the largest spin sector (with equal num- 
ber of spin-ups and spin-downs) is ( 18 ) = 48, 620. The 
large memory requirements for the calculation of the full 
spectrum, needed for the finite temperature NLCE, is 
the one that limits our expansion to the 9th order in this 
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TABLE I. Number of topologically distinct clusters up to the 13 th order of the expansion. 



order number 


1 


2 


3 4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


# of clusters 


1 


1 


2 3 


7 


14 


25 


50 


102 


216 


442 


929 


1969 



case. 

For the ground-state calculations, wc consider clusters 
up to the 13th order. Here, only in the last order, 1969 
clusters have to be diagonalized. Each of them has 26 
sites and the size of the Hilbert space for the largest spin 
sector is ( 13 ) = 10,400,600. The memory requirements 
for these calculations are much smaller than those for 
the full diagonalization (only three vectors of this size 
need to be stored at any given time). What limits our 



ground-state expansion to the 13th order is the time they 
take to converge in the Lanczos loop, so that the ground- 
state energy is obtained with a relative accuracy of 10~ 14 
for orders less than or equal to 12 and of 10 -9 for the 
13th (last) order. The convergence requirements that 
are imposed guarantee at least 5-digit accuracy for the 
final observables after summing the contributions of all 
of the clusters. 
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